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We present a method for optimizing quantum circuits architecture. The method is based on the 
notion of quantum comb, which describes a circuit board in which one can insert variable subcircuits. 
The method allows one to efficiently address novel kinds of quantum information processing tasks, 
such as storing-retrieving, and cloning of channels. 
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Quantum Mechanics plays a crucial role in the tech- 
nology of high precisions and high sensitivities, e. g. in 
frequency standards JDJ , quantum lithography H, two- 
photon microscopy |3y, clock synchronization [|4|, and 
reference- frame transfer jq| . In all these applications the 
essential problem is to achieve very high precision in: i) 
determining parameters; ii) executing a transformation 
that depends on unknown parameters. Since the param- 
eters are generally encoded with a transformation Q , as 
in the whole class of quantum metrology problems 0, 
and since the estimation itself can be considered as a 
special case of transformation (with classical c-number 
output), both tasks i) and ii) can be reduced to the gen- 
eral problem of executing a desired transformation de- 
pending on an unknown transformation. Taking into ac- 
count the possibility of exploiting N uses of the unknown 
transformation, the problem becomes to build a quantum 
circuit that has N circuits as input, and achieves the de- 
sired transformation as an output. This is what we call 
a quantum circuit board. 

A quantum circuit board is a network of gates in which 
there are N slots with open ports for the insertion of N 
variable sub-circuits (see Fig. Q). Since generally it is im- 




FIG. 1: A quantum circuit board. 

possible even in principle to achieve the desired transfor- 
mation exactly, the main task here is to optimize the cir- 
cuit board according to a given figure of merit — e. g. the 
gate-fidelity with the target transformation. A typical 
example is the optimal cloning of an undisclosed trans- 
formation U , which will be operated by a board with 
N slotted uses of U, and achieving overall in-out uni- 
tary transformation which is the closest possible to U® M 
with M > N. We emphasize that generally the over- 
all in-out transformation of the board and the slotted 
ones can be of any kind, including measurements and 
state-preparations (i. e. transformations to and from c- 
numbers), and that in some specific situations the slotted 
transformations can be even different from each other. 
In previous literature some special cases of circuit- 



board optimization have been considered, e. g. in re- 
gards to the discrimination of two unitary transforma- 
tions with N uses in parallel S, or in sequence M, and 
to the problem of phase estimation |Q, where an opti- 
mal quantum circuit architecture has been determined 
|l0| . However, the optimal architecture for estimating 
arbitrary unitaries is still unknown, and no systematic 
method for the general problem of circuit board opti- 
mization is available yet. 

In this Letter we present a complete method for opti- 
mizing the architecture of quantum circuit boards. After 
providing a convenient description of circuit connectiv- 
ity, we introduce the notion of quantum comb, which de- 
scribes all possible transformations operated by a quan- 
tum circuits board, and generalizes the notion of quan- 
tum channel to the case where the inputs are quantum 
circuits, rather than quantum states. We then present 
the optimization method, based on the convex structure 
of the set of quantum combs. The method allows one 
to reduce the apparently untractable problem of optimal 
circuit architecture to the optimization of a single pos- 
itive operator with linear constraints. Finally, we give 
two new applications in which the present approach dra- 
matically simplifies the solution of the problem. 

A quantum circuit operates a transformation from in- 
put to output, and is graphically represented by a box 
with input and output wires symbolizing the respective 
quantum systems. The quantum systems of different 
wires are generally different, and may also vary from 
input to output. Let us associate Hilbert spaces Hj n 
(H ou t) to all input (output) wires, and denote by p ul 
(Pout) the corresponding states. The action of the circuit 
is generally probabilistic, i. e. different in-out transfor- 
mations can randomly occur, as in a measurement pro- 
cess. Each transformation is described by a linear map 
Pin — * ^(pin) = kp out , with the proportionality factor 
k = Tr^p^)] sC 1 giving the probability that c € 
occurs on state p m . To describe a legitimate quantum 
transformation, the map H : Lin (Hj n ) — > Lin(H out ) Jll[] 
has to be completely positive (CP) Jiaj and trace non- 
increasing. Trace-preserving maps — i.e. deterministic 
transformations — are called quantum channels. Notice 
that a map ^ , rather than representing a specific cir- 
cuit, is univocally associated to the equivalence class of 
all circuits performing the same in-out transformation. 

The linear map c € can be conveniently rewritten using 
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the so-called " Choi- Jamiolkowski" representation p3| , 
corresponding to the following one-to-one correspondence 
between linear maps ^ : Lin(H; n ) — > Lin(H out ) and linear 
operators C € Lin(H out ® H in ) given by 



C = Choi(^) :=^®«/(|Q)(0|). 
%f(p) = Chor^CXp) :=Tr in [(/ out <: 



(1) 
(2) 



where is the identity map, \Q) is the unnormalized 
maximally entangled state = ^2 n \n)\n) 6 H® 2 , and 
T denotes transposition with respect to the orthonormal 
basis {\n)} for Hj n . The map ^ is CP if and only if the 
operator C — called Choi operator — is positive 14 . 



Two quantum circuits can be connected in all the ways 
allowed by the physical matchings between input and 
output wires (see e. g. Fig. ||, where the wires labelled 
d are connected): a connection will result in the compo- 
sition of the corresponding CP maps, and hence of the 
corresponding Choi operators. Since building a quantum 
network means connecting many circuits, it is crucial to 
have a handy way to describe circuit connectivity with 
minimum overhead of notation. We provide here three 
simple rules that accomplish this goal: 
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W = H * H 
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A = Choi(.c/), B = Choi(^), C = Choi(^) = A * B 

FIG. 2: Connection of two quantum circuits s4 and £3. Wires 
are labelled according to Rule |l| The Choi operator the re- 
sulting circuit H is given by the link product of Rule |^. 



Rule 1 (Labelling) Each quantum wire is marked with 
a different label, except for wires that are connected, 
which are identified with the same label. 

Rule 2 (Multiplication) The multiplication of two 
Choi operators A 6 Lin(H a ,k !Cf i) and B 6 Lin(Hd, e ./, 9 ) 

is intended in the tensor fashion, i.e. AB = {A ® 

Ie,f,g)(I a ,b,c®B). 

Rule 3 (Composition) The connection of two circuits 
with Choi operators A and B — acting on Hilbert spaces 
labelled according to Rule |J — yields a new circuit with 
Choi operator C given by the link product 



C = A*B = Trj[A 0J B], 



(3) 



6j denoting partial transposition over the Hilbert space J 
of the connected wires, and the multiplication in square 
brackets following Rule |^. 



Rule [| directly follows from Eqs. ([!]) and (||). Notice that 
due to invariance of trace under cyclic permutations, the 
link product is commutative: A * B = B * A. Using it, 
the action of a linear map ^ on a state p in Eq. (j|) 
can be rewritten as ^(p) = C * p. Assembling many 
circuits , ^2 , ■ • • , % yields a quantum network whose 
Choi operator is simply given by C = C\ * C2 * • • • * Ck ■ 
We are now ready to treat quantum circuit boards. To 
start with, we consider the case of a deterministic circuit 
board, i.e. a network of quantum channels with TV open 
slots for the insertion of variable subcircuits. It is clear 
that by reshuffling and stretching the internal wires any 
circuit board can be reshaped in the form of a " comb" , 
with an ordered sequence of slots, each between two suc- 
cessive teeth, as in Fig. ||. The order of the slots is the 
causal order induced by the flow of quantum informa- 
tion in the circuit board. We label the input systems 
(entering the board) with even numbers 2n, and the cor- 
responding output systems (exiting the board) with odd 
numbers 2n + 1, with n ranging from to N . 
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FIG. 3: Every circuit-board can be reshaped in form of a 
"comb", with an ordered sequence of slots, each between two 
successive teeth. The pins represent quantum systems, enter- 
ing or exiting from the board (the horizontal arrow represents 
the quantum information flow). 

A quantum comb with N slots is clearly equivalent to 
a concatenation of JV + 1 channels with memory, which is 
in turn equivalent to causal network, namely a network 
where the quantum state of the output systems up to 
time n does not depend on the state of the input systems 
at later times n' > n, with n, n' € {0, 1, . . . , N} The 
causal network can be easily obtained by redrawing the 
comb as an equivalent circuit with all inputs on the left 
and all outputs on the right, as in Fig. ||. We define the 
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FIG. 4: Each quantum comb is equivalent to a causal net- 
work, with all inputs on the left and all outputs on the right. 
The Choi operator of a comb is the Choi operator of the cor- 
responding causal network. 
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Choi operator of a quantum comb as the Choi operator 
R of the corresponding causal network. In terms of the 
Choi operator R, causality is equivalent to a set of linear 
constraints 



Tr 2 „ 



+ 1 



R {n) = hn®R {n ~ 1] , n = 0, 
RW = R, R^ = l, 



■ N, 



(4) 



where Tr2 n +i denotes the partial trace over the Hilbert 
space Hbn+i of the output wire labeled 2n + 1, Im the 
identity operator over the Hilbert space H2 n of the input 
wire labeled 2n, R {n) = Choi(<*f( n )), and 'afM i s the map 
of the (n + 1)- subnetwork from the first n + 1 inputs to 
the first n + 1 outputs. Precisely, we have the following 

Theorem 1 Every positive operator ^ i? G 
Uri(<8^=o~ 1 H 7 -) satisfying the linear constraints fy), is the 
Choi operator of a deterministic quantum comb. 

Proof By definition, it is enough to show that any 
operator R ^ normalized as in Eq. (^) is the 
Choi operator of a causal network. A causal network 
with N + 1 input/output pairs is described by a fam- 
ily of channels < «f( n ', n — 0, 1, . . . , N with the property 
Tr 2n+1 [<«f (")(/>("))] = tfC"- 1 ) (Tr 2 „[p(™)]), for any state 
p( n ) of the first n + 1 input systems. Using the corre- 
spondence of Eq. (|2|), one can easily see that this is 
equivalent to the normalization of Eq. (^) . 

A quantum comb transforms a series of N input cir- 
cuits , . . . , ^jv into an output circuit c €' depending on 
them (Fig. ||a). This transformation of circuits cor- 
responds to an TV-linear CP-map that sends the input 
Choi operators into the ouput Choi operator according 
to C = C\ * ■ ■ ■ * Cn * R, with R the Choi opera- 
tor of the comb. We call the mapping between circuits 
{^i, . . . , ^jy} !— > c &' supermap as it sends channels into 
channels, rather than states into states. Notice that, 
depending on the number of slots that are saturated a 
quantum comb can transform a series of circuits into a 
comb (Fig. ||b), or, more generally, a comb into a comb 
(Fig. ||c). As a matter of fact, a quantum comb realizes 
many possible mappings, all obtained by the link prod- 
uct with its Choi operator R. Therefore, the quantum 
comb can be completely identified with its Choi opera- 
tor. Remarkably, also the converse is true: any abstract 
supermap sending channels into channels in a CP fashion 
can be physically realized by a quantum comb Jl6| . 

The tools presented above provide a powerful method 
for optimizing quantum circuit architecture. Suppose we 
want to design a circuit board maximizing some convex 
figure of merit, e.g. the fidelity of the ouput circuit ^" 
with a desired unitary gate °U '. In our framework the 
optimization of the board architecture is reduced to the 
search of the optimal operator R ^ with the linear 
constraints (Q) . This is a standard problem of convex op- 
timization, for which efficient algorithms are known. Ba- 
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FIG. 5: A quantum comb realizes different transformations of 
quantum circuits, namely it can send (a) A series of channels 
into a channel, (b) A series of channels into a comb, or (c) 
An input comb into an output comb. 



sically, we only need to implement the search on the ex- 
tremal points of the convex set of Choi operators. More- 
over, the complexity of the search can be dramatically 
reduced by exploiting additional constraints, e.g. sym- 
metry properties of the circuit board. The optimal Choi 
operator will finally single out the optimal architecture, 
automatically deciding if the N slots of the circuit board 
have to be connected in a causal order or in parallel, or 
in any combination of the two. 

We illustrate our method in two concrete applications. 
The first application is the optimal universal cloning 
of unitary transformations, i. e. the problem of de- 
signing a quantum board that optimally achieves the 
N — > M cloning of an unknown unitary U G SU(d) 
in dimension d. The board has N slots containing N 
identical uses of the unknown unitary U and performs 
a transformation which is the closest possible to JJ® M . 
Using channel fidelity as the figure of merit, the prob- 




-m- 



FIG. 6: Cloning of unitary transformations. A quantum 
board with N = 3 input slots, designed to optimally emu- 
late M = 5 uses of the unknown unitary transformation U 
with maximum channel fidelity. 

lem is to find the Choi operator R that maximizes the 
average over all unitaries of the overlap Fjj (N, M) = 
l/{d 2M+N ){U\® M {U*\® N R \U)® M \U*)® N , where \U) = 
(U ® with = Y, n \n)\n). The architecture op- 

timization is then reduced to a standard convex analy- 
sis problem. For N = 1 and M = 2, we derived the 
optimal quantum board, achieving fidelity F "(1,2) = 
(d + ^/d 2 — l)/d 3 , significantly higher than the classical 
threshold reached by the optimal estimation of a unitary 
F est (l,2) = 6/d 4 for d > 2, F esi (l,2) = 5/16 for qubits 
|17| , this showing the advantage of coherent quantum in- 
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formation processing over any classical cloning strategy. 

Another interesting application is the storage and re- 
trieval of an undisclosed unitary transformation U from 
N uses, also called optimal quantum- algorithm learning. 
The problem arises from the need of running an undis- 
closed algorithm (available for N uses) on an input state 
tjj which will be available at later time. To this purpose 
one can slot the N uses of U in a quantum circuit board, 
put the output state of the board in a quantum memory, 
and, when the input state will be available, use the mem- 
ory to recover the unitary. The series storing-retrieving 
is represented by the quantum comb in Fig. ^ which can 
be cut into two parts, a storing one including only the 
uses of U, and a retrieving one including ip (the output 
state of the first part is stored in a quantum memory 
and is then fed in the second part). Also in this case 
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Fig. Again, by optimizing the operators {Ri} one 
will automatically determine the optimal disposition of 
the unitaries in the circuit, a problem whose solution is 
up to now known only in the very special case of phase 
estimation [10|. 



FIG. 7: Quantum-algorithm learning. One wants to run an 
undisclosed unitary U on a quantum state t[> which is available 
after the lapse of time in which the uses of U are available. 

our method reduces the optimization to a convex anal- 
ysis problem. For N = 1 and N = 2 we found average 
gate-fidelities F = -jp and F = -p, respectively, which 
coincide with the value attained by the optimal estima- 
tion of unitaries |l§f| . Remarkably, the optimal universal 
storing/retrieving of unitaries does not need a coherent 
interaction with the quantum memory at the retrieving 
stage — which is purely classical — , but only an entan- 
gled input state at the storing stage. This means that 
the quantum memory is not really needed, and that the 
learned algorithm can be executed an unlimited num- 
ber of times with constant performance. The situation 
is radically different if no entanglement is allowed in the 
storing stage: in this case the optimal retrieving is purely 
quantum, yielding an optimal learning that is forgetful. 

We conclude by mentioning the extension of our 
method to the optimization of probabilistic circuit 
boards, containing measuring devices that produce dif- 
ferent transformations depending on random outcomes. 
The probabilistic comb corresponding to outcome i will 
have Choi operator Ri, with the sum over all outcomes 
Y], Ri = R giving the Choi operator of a deterministic 
comb. Indeed, introducing a classical register with or- 
thogonal states \i) € He we can define R = Y2i Ri ® 
which is the Choi operator of a deterministic comb with 
H2W+1 := H2Ar+i <8> He- The comb corresponding to Ri 
is then obtained after applying the comb of R, by mea- 
suring the register on the basis {\i}} and postselecting 
outcome i. Probabilistic combs are a fundamental tool 
to address the optimized circuit architecture for estimat- 
ing unknown transformations with multiple copies (see 




FIG. 8: Quantum comb for the estimation of unitary trans- 
formations with multiple uses. 
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